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A mathematical model describing the coupling between two independent amplification mecha- 
nisms in auditory hair cells is proposed and analyzed. Hair cells are cells in the inner ear responsible 
for translating sound-induced mechanical stimuli into an electrical signal that can then be recorded 
by the auditory nerve. In nonmammals, two separate mechanisms have been postulated to contribute 
to the amplification and tuning properties of the hair cells. Models of each of these mechanisms 
have been shown to be poised near a Hopf bifurcation. Through a weakly nonlinear analysis that 
assumes weak periodic forcing, weak damping, and weak coupling, the physiologically-based models 
of the two mechanisms are reduced to a system of two coupled amplitude equations describing the 
resonant response. The predictions that follow from an analysis of the reduced equations, as well as 
performance benefits due to the coupling of the two mechanisms, are discussed and compared with 
published experimental auditory nerve data. 
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I. INTRODUCTION 



The natural environment presents the auditory system with the challenge of responding to sounds over many orders 
of magnitude; the threshold of hearing and the threshold of pain differ by about thirteen orders of magnitude. For 
the ear to discriminate between sounds over such a large dynamic range, it is necessary for auditory stimuli to be 
compressed into a much smaller, more achievable range of physical responses. This is accomplished through a nonlinear 
mechanism in which small amplitude sounds are amplified to a greater extent than larger amplitude sounds. In order 
to process complex sound stimuli, it is also necessary for the auditory system to distinguish between the frequency 
components of the stimuli. The ear's amplification and frequency discrimination properties are thought to be derived 
. from a common, metabolically-powered mechanism, the details of which have been the topic of much investigation 

o: ELII- 

In both mammals and nonmammals, hair cells of the inner ear are responsible for translating sound- induced 
mechanical stimuli into a neurotransmitter signal which induces the firing of the auditory nerve [j, Q. Each hair 
cell consists of a cell body which is contacted from below by the auditory nerve and a hair bundle consisting of 
actin-supported fibers. When sound stimulates the auditory organ, the resulting motion of the hair bundle causes 
transduction channels to be mechanically pulled open. Ionic current then enters the cell body through the transduction 
channels, thereby depolarizing the cell, and ultimately causing the release of neurotransmitter at the auditory nerve 
synapse. 

In mammals, the frequency-discrimination properties of the basilar membrane, the membrane in which the hair cells 
are embedded, contribute to the auditory system's capacity to distinguish between sounds of different frequencies. By 
contrast in nonmammals, the surface in which the hair cells are embedded lacks tuning properties. The nonmammalian 
auditory system is thought to achieve its frequency tuning properties through two different mechanisms, both intrinsic 
to the hair cell. The first mechanism involves the mechanical motion of the hair bundle. Experiments indicate that 
the hair bundle responds actively, with greater energy than provided by the stimulus, if forced near its resonance 
frequency Q. Evidence for a second mechanism, referred to as electrical resonance, is provided by the decaying 
oscillations that are observed in the membrane potential of the cell body in response to constant current injection 0]. 
These oscillations indicate that the cell body possesses a preferred response frequency. 

Dynamical systems methods have proven useful in analyzing the frequency tuning and amplification properties of 
physiologically-based auditory models. Interestingly, models of both the active motion of the hair bundle and the 
electrical resonance mechanism have been shown to be poised near a Hopf bifurcation 0, @. A Hopf bifurcation is 
a robust mechanism for generating spontaneous oscillations as a control parameter of a nonlinear system is varied. 
It occurs when a static equilibrium loses stability via a complex conjugate pair of eigenvalues (of the associated 
linear stability problem) crossing the imaginary axis in the complex plane with nonzero imaginary part. It has been 
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suggested that the hair cell critically tunes itself so that its parameters are poised just below the bifurcation point, 
thereby making the cell sensitive to stimuli at the Hopf bifurcation frequency, without causing spontaneous oscillations 
[ToL [Hi . These investigations determine the generic frequency tuning and amplification properties of a periodically- 
forced system in the vicinity of a Hopf bifurcation; specifically, they analyze the characteristics of solutions that are 
frequency-locked to a weak, additive resonant forcing term. Sufficiently close to the Hopf bifurcation point, the system 
is compressively nonlinear: small inputs are amplified to a greater extent than larger ones [Tel lllj. Moreover, the 
compression of the dynamic range is accompanied by frequency tuning, which is sharper for small amplitude inputs 
than for larger amplitude signals. 

Previous studies of hair cell amplification models have considered the normal form for a system near a Hopf 
bifurcation without actually performing the normal form reduction from the physiologically relevant mathematical 
model. Here, in Appendix II, we reduce the Hudspeth and Lewis model of the electrical tuning mechanism [12, EH 
to the normal form for a system near a Hopf bifurcation, thereby determining the numerical values of the coefficients 
of the normal form corresponding to the model and parameters used by Hudspeth and Lewis. We find that the 
coefficient of the nonlinear term in the normal form has comparable real and imaginary parts; it is not purely real as 
was assumed in earlier investigations [IcL . We show that a result of the nonzero imaginary part is that the response 
of the system to resonant forcing may be hysteretic and that the frequency tuning curves are no longer symmetric 
about the resonance frequency. We further propose a model that describes weak coupling between the hair bundle 
amplification mechanism and the electrical resonance mechanism. We assume that both oscillation mechanisms are 
critically tuned to approximately the same resonant frequency and only weakly damped so that the Hopf bifurcation 
normal form applies to each independently, i. e. when the other mechanism is suppressed. We then assume weak linear 
coupling of the mechanisms, and direct forcing of the hair bundle at a frequency that is close to its natural frequency. 
As in earlier investigations, the analysis focuses on the frequency-locked solutions, and, in particular, on how the 
magnitude of response grows with the forcing. We find that the combined critically-tuned amplification system, can 
lead to a response, R, that scales with F 1 / 9 , where F is the resonant forcing amplitude, thereby leading to enhanced 
amplification, R/F, of small signals. We also explore the enhanced frequency-tuning characteristics of the combined 
mechanical and electrical amplification system, comparing it with those associated with a single tuning mechanism. 

Our paper is organized as follows. In section II we introduce the reduced mathematical model, with the mathematical 
details of the reduction from the physiologically-detailed models relegated to Appendix II. Section III contains our 
analysis of the reduced model, focusing particularly on the simpler situation of unidirectional coupling from the hair 
bundle tuning mechanism to the electrical resonator. We present response- versus-forcing curves that demonstrate the 
transition from a linear response (R oc F) to a response R cx F 1 / 9 as the amplitude of the (weak) signal increases. 
We also demonstrate the sharper tuning that is possible with the combined amplification system. Finally, section IV 
compares our model predictions with published experimental auditory nerve data. 

II. MODEL 

Two distinct mechanisms contribute to auditory tuning in nonmammalian vertebrates: an 'electrical' resonance 
arising from an interplay between ionic currents through the cell membrane, and a 'mechanical' resonance associated 
with the active motion of the stereocilia in reponse to stimuli at their resonance frequency. We start our discussion 
by focusing on the electrical resonance mechanism. The underlying biophysical components have been discussed by 
Hudspeth and Lewis [l2l . [l3| , who performed a set of experiments that carefully characterized the dynamical properties 
of the major ion channels on the cell bodies of bullfrog saccular hair cells. On the basis of these experiments, they 
developed a single compartment model of the hair cell (see Appendix I) using the simplifying assumption that only 
two major active ion channels, a voltage-gated calcium channel and a calcium-gated potassium channel, contribute 
to the cell's dynamical behavior. In this model, the dynamical evolution of the membrane potential V m is given by 
an equation based on the direct application of Kirchoff 's laws to a circuit that represents the flow of ions across the 
membrane: 



- C m ^- = 9Cam 3 (V m - E Ca ) + 5K(Ca)(0 2 + 3 )(V m - E K ) + g L (V m - E L ) - I. (1) 

Here V m is the membrane potential and C m is the membrane capacitance per unit area. The voltage-gated calcium 
(Ca) current is represented by g caXn (V m — £?Ca)) where gca is the maximum Ca conductance per unit area, m is 
the voltage-dependent fraction of open conformational subunits in the Ca channels, and i?ca is the reversal potential 
for the Ca ion channels. The Ca-gated potassium (K) current is represented by gK(Ca){02 + 3 ){y m — E%), where 
<?K(Ca) is the maximum K conductance per unit area, (O2 + O3) is the fraction of K channels in one of their two 
open states, and Ek is the reversal potential for the K ion channels. The term gL(V m — El) represents all passive 
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FIG. 1: Numerical simulations of the response of the membrane potential of a hair cell in the Hudspeth and Lewis model [12l.[l3|. 
to a constant current, /, injected at t — .06. (a) I = 65 pA; (b) I = 95 pA. A Hopf bifurcation occurs at I* « 91.3 pA. The 
other parameters of the H&L model (|15p . used in the simulations, are given in Appendix I. 



ion channels as a leak conductance ql per unit area and a reversal potential El- The command current / is directly 
injected into the cell body. The formulation of the model involves six additional equations (see Appendix I) that 
describe the dynamical evolution of the fraction m of open units in the Ca channels, the intracellular concentration 
of Ca ions close to the cell membrane, and the fraction of Ca-gated K channels in each of their three closed states 
(Co, Ci, C2) and two open states (O2, O3). In their seminal work [12, EH, Hudspeth and Lewis (H&L) experimentally 
characterized the value of the various parameters that appear in these equations. 

The H&L model reproduces qualitatively the decaying membrane potential oscillations observed in current-clamp 
experiments in which a current of constant amplitude / is injected into the cell body [HI, [l3j]. As shown in Figure [TJ 
the response of the membrane potential to a step current of amplitude / depends crucially on the value of the injected 
current. For small injected currents, as illustrated in Figure[TJa), the membrane potential exhibits an oscillatory decay 
to a new constant value. For larger current values, as illustrated in Figure []Jb), the membrane potential decays to a 
new state that is oscillatory. This qualitative difference signals a transition between a regime in which the asymptotic 
state is a fixed point and a regime in which the asymptotic state is a limit cycle. In the H&L model, this transition 
occurs by a Hopf bifurcation [9j . 

A Hopf bifurcation occurs when a fixed point of a system of ODE's undergoes a change in stability in which a 
complex conjugate pair of eigenvalues A, A passes from the i?e(A) < to the i?e(A) > side of the imaginary axis in 
the complex plane. Figure [2] shows the evolution of the eigenvalues of the H&L model linearized around its fixed point, 
as the input current is increased from and 100 pA. Note that three of the eigenvalues are always real and negative, 
and one complex conjugate pair remains in the Re(X) < semispace. The leading complex conjugate pair crosses the 
Re(X) = axis for I* ss 91.3 pA; this is the value of the input current at which the fixed point becomes unstable, 
as determined by fixing the parameters of the model at the experimentally-based estimates listed in Hudspeth and 
Lewis's original paper. That the H&L model is, for physiologically reasonable values of the parameters, poised near a 
Hopf bifurcation has profound implications for the signal processing capabilities of the modelled hair cell. Specifically, 
at the Hopf bifurcation the system is compressively nonlinear, as it exhibits a large amplification of small amplitude 
inputs and a smaller amplification of large amplitude inputs [T5L [TTj] . Moreover, the resulting compression of the 
dynamic range is accompanied by a sharp frequency tuning of the response to small amplitude inputs and a broad 
tuning in the response to large amplitude inputs [lfj, [ll| . 

The dynamical behavior of the model is asymptotically described by the amplitude of the mode associated with the 
most unstable eigenvector, which is the one associated with the complex conjugate pair of eigenvalues that cross the 
Re(X) — axis at the Hopf bifurcation. Sufficiently close to the bifurcation, the system can be generically described 
by a normal form equation [14| : 



— - (a + ib)A+(c + id)\A\ 2 A, (2) 

which characterizes the dynamical evolution of the complex amplitude A of the most unstable mode. In this equation, 
the parameter a is a measure of the distance to the Hopf bifurcation at a = 0; a is negative below the bifurcation and 
positive above. The parameter b is the frequency of the system at the bifurcation, and the parameter d measures the 
shift in preferred frequency as the amplitude of the solution increases, as readily seen by setting A = re* n * in to 

establish that r = . and Q = b+ dr 2 . The parameter c distinguishes between supercritical (c < 0) and subcritical 



1500r 

1000- 
500- 

¥ 

-500- 
-1000- 
-150 



0000 -8000 -6000 -4000 -2000 2000 

Re(X) 



FIG. 2: The eigenvalues of the Hudspeth and Lewis model (|15|l . linearized about its static equilibrium state as described in 
Appendix II, evolve in the complex plane as / is increased from pA to 100 pA. Solid lines indicate the trajectory of the 
eigenvalues with increasing I. Other parameters of the equations are given in Appendix I. 

(c > 0) Hopf bifurcations. 

We have carried out a standard nonlinear reduction of the H&L model to the normal form (see Appendix II). 
In the resulting normal form ()26|) . a is proportional to AI = (7 — I*)/ 1* and c is negative. This reduction thus 
establishes the supercritical nature of the Hopf bifurcation in the H&L model. For this supercritical bifurcation, there 
is a transition from a stable fixed point to a state in which the fixed point becomes unstable and a stable limit cycle 
exists. Moreover, in the supercritical case, the radius r of the limit cycle grows, with distance past the bifurcation 
point, with the characteristic scaling r ot \J AI. 

In the current-clamp experiments conducted by Hudspeth and Lewis [ill IHj], the current injected into the cell was 
of constant amplitude. In contrast, natural inputs to the cell are due to a time-dependent, sound-induced mechanical 
displacement of the hair bundle, which results in time-dependent changes in the conductance through the transduction 
channels. Because the transduction channels are passive, this time-dependence can be incorporated into the H&L 
model through the leakage conductance term on the right-hand-side of (jT|) . In the case of a simple time-periodic 
conductance, the reduction carried out in Appendix II suggests that for command currents close to I* the most 
important contribution of the time-periodic forcing to the asymptotic dynamics comes from the Fourier component 
that is closest to the resonance frequency of the system. In the case of a weak periodic signal with frequency close 
to the resonator's frequency, the time-periodic input can be represented as an additive contribution to the amplitude 
equation, which then takes the form: 

dA 

— = {a + ib)A+{c + id)\A\ 2 A + Fe luJt . (3) 

Due to the time-translation symmetry of the unforced case, we can, without loss of generality, assume that F is real 
and positive. We note that earlier investigations of |[3J in the context of amplification mechanisms in auditory hair 
cells [ToL lllT | assumed a real coefficient of the nonlinear term, i. e. they made the nongeneric assumption that d = 
such that the preferred frequency of the nonlinear oscillator had no amplitude-dependence. 

We next consider the tuning mechanism associated with the mechanical deflection of the hair bundle. Experiments 
in which a glass fiber was attached to the hair bundle and used to mechanically stimulate the bundle at a specific 
frequency showed that the hair bundles respond preferentially to stimuli at their resonance frequency [6| . The motion 
of the hair bundles has been shown to be sensitive to the amount of calcium ion entering the transduction channels 
[l5| . A model for hair bundle motion due to calcium binding within the stereocilia was proposed by Choe et al. Q . In 
their model, when transduction current enters the hair bundle Ca ions attach themselves to the transduction channels 
at sites within the stereocilia; this attachment causes an increase in tension, which in turn causes the channel to close 
Q. When the transduction channel closes, the Ca ions detach themselves from the binding sites and the channel 
returns to its regular tension allowing the cycle to repeat itself. Choe et al. have shown that, when the parameters of 
their model are specified within physiologically reasonable ranges, the model is near a Hopf bifurcation. Experimental 
evidence also indicates that the relationship between stimulus magnitude and magnitude of the hair bundle oscillations 
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obeys the scaling that would be expected for a system tuned near a Hopf bifurcation [16(, lending support to our 
assumption that a model describing the hair bundle dynamics is poised near a Hopf bifurcation. Assuming this second 
mechanism is tuned sufficiently close to a Hopf bifurcation, it too can be reduced to the normal form {2J for a system 
near a Hopf bifurcation. 

It only remains to consider first the manner in which the two tuning mechanisms are coupled in the biological 
system, and second how this coupling should be represented in the reduced model. One clear source of coupling 
between the two tuning mechanisms is through the transduction current. The magnitude of the transduction current 
entering through the stereocilia is directly related to the magnitude of displacement of the stereocilia. This relationship 
between displacement and the amount of current entering the cell has been measured indirectly by measuring the 
change in the receptor potential of the cell in response to stereocilia displacements of different magnitudes [l7l ]. 
Such experiments indicate that, in absence of stimuli, a small amount of current is flowing into the cell through the 
stereocilia. When the hair bundle is deflected in the negative direction, away from the tallest stereocilia, transduction 
channels close and the amount of current flowing into the cell decreases. Similarly, when the hair bundle is deflected 
in the positive direction, the amount of current entering through the transduction channels increases and eventually 
saturates. For small displacements in the positive direction, the relationship between the displacement of the hair 
bundle and the change in the receptor potential is approximately linear [17| . As the amplitude of the hair bundle 
oscillations increases due to stimulation at the stereocilia's resonance frequency, the amount of current entering the 
cell body and providing a forcing to the second electrical resonance mechanism increases. This clearly provides a 
means of coupling from the stereocilia tuning mechanism to the electrical resonance mechanism. 

There is also evidence for coupling in the reverse direction, from the electrical resonance mechanism to the hair 
bundle resonance mechanism, in that electrical stimulation of the cell body has been shown to induce displacement 
of the stereocilia [l8l [T9L I20I [2ll |22| . The exact mechanism for coupling in this direction is less clear. Experimental 
comparisons between current injected into the cell body and the resulting displacement of the stereocilia indicate that 
the linear approximation is reasonable for small current injections [23]. The presence of coupling in both directions 
raises the question of whether there are actually two separate tuning mechanisms or the stereocilia tuning mechanism 
is simply a manifestation of the electrical tuning mechanism. This question was addressed in experiments in which 
the cell body was voltage clamped to silence electrical resonance oscillations allowing the motion of the stereocilia 
to be probed separately [24J. These experiments demonstrate active motion of the stereocilia even in the absence of 
the electrical resonance mechanism. Electrical resonance experiments are often performed by direct current injection, 
with transduction channels blocked, so the independence of the electrical resonance mechanism from the active motion 
of the stereocilia was never in question. 

From the biological evidence, it is reasonable to assume that the coupling between the two mechanisms is linear 
for sufficiently small forcing. This leads us to a reduced model consisting of two coupled amplitude equations of the 
form, 

— - = ( ai + ib x )Ai + (ci + id 1 )\A 1 \ 2 A 1 + 7l e*A 2 + Fe^\ (4) 



— - = (a 2 + ib 2 )A 2 + (c 2 + id 2 )\A 2 \ 2 A 2 + l2 e^A x . (5) 
at 

In this model, Q represents the hair bundle resonance mechanism which receives a sound-induced time-dependent 
forcing (cx F) as well as feedback from the electrical resonance mechanism (cx A 2 ). Equation (J5]) represents the 
electrical resonance mechanism which receives a forcing proportional to the displacement of the stereocilia. Note that 
we allowed for a phase difference ipj in each of the coupling terms, with the corresponding parameters 7,- taken to be 
real and non-negative. Moreover, we note that when both 71 and 72 are nonzero, then the Hopf bifurcations that cause 
spontaneous oscillations in the unforced problem (F — 0) will shift away from a,j = 0. As detailed mathematically 
in Appendix II, the model is valid for the case in which each system is tuned close to the Hopf bifurcation (|aj| 
sufficiently small), each system is tuned near the resonance frequency (bj ~ lu), and the forcing and the coupling are 
weak (F, 7 j sufficiently small). 

We have determined from Hudspeth and Lewis's model the numerical values of the coefficients a 2 , b 2 , c 2 , d 2 of 
for H&L's physiological parameters. Models also exist for the stereocilia mechanism, so it is possible that the same 
coefficients of Q could be determined based on these models. Performing the second reduction however would not be 
particularly useful in this paper because our analysis requires that both systems be tuned close to the same frequency. 
There are multiple ways to tune the physiological models such that they yield vibrations at a required frequency. 
Thus, without a very clear idea of the physiologically reasonable method to adjust the model parameters, there is 
no way to determine a consistent result for the numerical values of the coefficients in the amplitude equation ([4]). 
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Nonetheless, many of our conclusions regarding scaling laws hold provided that the Hopf bifurcations are supercritical 
(i.e. provided Cj < 0), and that our fundamental modeling assumptions are met. 



III. ANALYSIS 



A. Response- Versus-Forcing Relationship 

Our analysis of the reduced model, Eqs. ((I])-©, focuses on frequency-locked solutions of the form Aj = Rje^ t+ ^ j \ 
j = 1,2, where Rj > and 4>j G [0, 27r) are constants. We wish to determine how the magnitude of the electrical 
response, measured by R 2 , scales with the sound-induced mechanical forcing amplitude F. This scaling depends on 
the proximity to the Hopf bifurcation, captured by the linear damping coefficients aj < in ©-©J. It also depends 
on how closely tuned the natural frequencies, bj, of the nonlinear oscillation mechanisms are to each other and to the 
driving frequency, u>. 

Substituting Aj = R J e l( ' UJt+ ^'> into Q-© yields the following pair of complex-valued algebraic equations defining 
an implicit relationship between the real quantities F and R 2 : 

Fe-^ 1 = -{ai+i{b 1 -oj))R 1 -{ci+id 1 )Rl~ 1 iR 2 e i{ ' t>2 -' t,1+x l' l \ (6) 
72 iJ 1 e < {* 1 -* a+ *^ = -{a 2 +i(b 2 -u))R 2 - (c 2 +id 2 )R 3 2 . (7) 



We can solve © for i?i in terms of R 2 : 

e — 1(^1-^2+^2) 
i?i = 



?2 



((a 2 + t(b 2 - uj))R 2 + (c 2 + id 2 )R$\ . (8) 



Here we can determine the phase difference (4>i—(t> 2 ) by the requirement that R\ be real and nonnegative. Substituting 
this expression into ([6]), we find 

F = e 1 ^ 2 -^ f ai R 2 + a 3 Rl + a 5 R% + a 7 R 7 2 + a 9 R^ , (9) 

where 

a x = - lie ^+^ + — ( ai +i(h-oj))(a 2 +i(b 2 -u)), 



a 3 



72 

— (ca + id 2 )(ai + i(h - «)) + i- e - 2l (^-^+^)( Cl + id 1 )(a 2 + z(6 2 - c^)) 3 , 



72 ' 72 



, lr _ = Ae- 2 ^ 1 - ^,2+ ^)(a2+^(62-^)) 2 (cl+^dl)(c2 + ^ d2), (10) 



7 2 3 
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7 2 

J_ e -2 4 (0 1 -0 2+ V 2 ) (ci + idl)(ca + id2 )3 _ 
72 



Again, the phase </>2 in ([9]) is determined by the requirement that the forcing magnitude F be real and non-negative. 
It immediately follows from the polynomial form of @ that the response R 2 need not be a single-valued function of 
F . It also follows that, with increasing forcing, there is a transition from a linear regime, R 2 cx F, for sufficiently 
small forcing, to a regime where R 2 oc i^ 1 / 9 for larger values of F. However, whether this transition occurs for small 
values of F, for which the model (JU)-© is valid, depends on both the magnitude of the coupling coefficients 7^ and 
the magnitude of the linear coefficients aj + i(bj — oj). In particular, if we let e — \a 2 + i(b 2 — uj)\, then we expect the 
transition to the regime R 2 cx F 1 / 9 will occur for small F provided that |ai + i(b\ — w)|7| is at most 0(e 3 ), and 717! 
is at most 0(e 4 ). 

The response vs. forcing characteristics associated with the model (@|-([5]) are further explored in Fig. [3l For this 
figure, we assume supercritical Hopf bifurcations associated with the mechanical and electrical resonance mechanisms 
so that Cj < in the model equations. With an appropriate re-scaling of amplitudes Aj, we may then assume Cj = — 1 
for both j = 1 and j — 2. Our direct calculation of the normal form coefficients in ([S]) (see Appendix II), from the H&L 
model, yields d 2 /c 2 ~ 1.1. Thus in many of our numerical computations we set d 2 = —1.1. Finally, if we scale time by 
the dimensioned forcing frequency, we may take u> = 1 in the model equations. Fig. [^indicates stable frequency-locked 
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solutions as solid lines and unstable solutions by dotted lines. The linear stability of the frequency-locked solutions is 
determined by substituting the ansatz Aj = Rje i ^ u ' t+ ^ (1 + Zj(t)) into Eqs. (HI)-© and then linearizing about Zj = 0. 
We then find that the perturbations Zj satisfy the following system of linear differential equations 
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where z\ denotes the complex conjugate of z\, etc., and 



Mi 
M 2 

M 3 



a± + ib\ — ilu + 2(ci + id±) 
(ci + ie?i)i? 2 



7i 



^ 2 r H4>i-4>i ) 



M 4 = a 2 + ib 2 - iuj + 2(c 2 + id 2 ) 
M 5 = (c 2 +id 2 )R 2 2 

Me 



(12) 



#1 '((01-02+^2) 



The solution Aj = Rje lt - ut+t ^ j \ with i?j and <f>j satisfying ©-fT]), is stable if the eigenvalues of the matrix associated 
with the linearized problem all have negative real part. 

Figure [3]^ a) demonstrates the predicted transition from linear response to nonlinear response with R 2 oc F 1 / 9 . 
Figure [3][b) shows that as the damping in the system increases, larger forcings are necessary to reach this nonlinear 
regime. Moreover, this plot demonstrates the possibility of hysteresis which can occur when bj — u> and dj have 
opposite signs, provided the damping is not too great. Figure [3fc) shows how the transition from linear to 1/9 scaling 
moves to larger forcing when the detuning b 2 — ui is increased; note that in this plot a scaling of R 2 cx F 1 / 3 is also 
evident over an intermediate range of forcings. Changes in the magnitude of the feedback coefficient 71 can either 
enhance or degrade the response, R 2 , depending upon the coupling phases ipi and ip 2 . Figure G2d) shows an example 
of the change in the response versus forcing relationship as the feedback coefficient, 71, is increased. 



B. Uniqueness and Stability: Unidirectional Coupling 

The problem of determining the uniqueness and stability of solutions is greatly simplified in the case of unidirectional 
coupling between the mechanical and electrical resonators. In this section we focus our further analysis on the case 
where the feedback from the electrical resonator to the stereocilia can be neglected, i.e., we focus on the case 71 = 
in gj). 

In this unidirectional-coupling case, M 3 = in the stability matrix (|11[) . and the linear stability problem simplifies 
to one of determining the eigenvalues associated with the 2 x 2-blocks on the diagonal. For instance, if we let <j\ and 
02 be the eigenvalues associated with the [z\, zi)-equations, we find 

01 + cr 2 = Mi + Mi = 2(oi + 2ciR\) (13) 
a X G 2 = \Ah\ 2 -\M 2 \ 2 = aj + {bi - uj) 2 + 4[a lCl + {h - uj^Rj + 3(4 + dj)Rf. 

Similar equations for the remaining two eigenvalues, associated with the (z 2 , "^-equations, hold: they are obtained 
from (Tl3|) by interchanging the 1 and 2 subscripts on its right-hand-side. In the case of supercritical Hopf bifurcations 
(cj < 0), and damping of spontaneous oscillations (aj < 0), the frequency- locked solutions are stable provided 
|Mi| 2 — |M 2 | 2 > and |M 4 | 2 — |M 5 | 2 > 0. While these conditions hold for sufficiently small and sufficiently large 
amplitudes R\ and R 2 , they may be violated in the intermediate regime if a,jCj + (bj — uj)dj < for either j = 1 or 
j = 2. Since we are interested in the case that aj and Cj are both negative, a necessary condition for a frequency-locked 
solution to lose stability, with increase in forcing F, is that bj — to and dj must have opposite signs. To gain some 
insight into this criterion, it is useful to note that the preferred frequencies, 61 + d\R\ and b 2 + d 2 R\, of each tuning 
mechanism, in absence of forcing, are dependent upon response magnitude. If the natural frequency of the cell shifts 
away from the forcing frequency with increasing response ([bj — cu]dj > 0), then ajCj + (bj —u>)dj > and the entrained 
solution is always stable and unique. Instabilities, and their associated hysteresis in the reponse vs. forcing curves, 
can only occur in the unidirectionally-coupled case if the preferred frequency of at least one of the tuning mechanisms 
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FIG. 3: Sample log-log plots of response (R 2 ) versus forcing (F) from ([9l ll0[) . for the system with Ci = C2 = —1 and 
with varying damping, coupling, and detuning magnitudes. Solid (dotted) lines represent (un)stable frequency-locked so- 
lutions, (a) Transition from linear response (R 2 oc F) to nonlinear response with R oc F 1 ^ 9 ; dashed lines, for compar- 
ison, have slopes 1 and 1/9 as indicated. Parameters set to a\ = —0.02, a 2 — —0.001, 61 — ui = 0, b 2 — w = .01, 
d\ — —2, d 2 — —1.1, 71 = 0, 72 = .1, and ip2 = 0.647T. (b) Response curves obtained with linear damping parameters 
(ai,a 2 ) = (-.0002, -.0001), (-.0002, -.1), (-.2, -.1); other parameters set at 61 - u = .01, 62 - w = .02, d a = -6, d 2 = -4, 
71 = 0, 72 = .1, and ip2 — 0.647T. (c) Response curves obtained with detuning b 2 — uj = .001, .01, .1; other parameters set at 
61 — u = .002, ai = —.001, a 2 = —.002, di = —2, d 2 = —1.1, 71 = 0, 72 = .1, and ip 2 = 0.647T. (d) Response curves obtained 
with different backward coupling magnitudes of 71 = 0, .1, .2, .3, .4; other parameters set at a\ = —0.1 a 2 = —0.2, 61 — us — 0.02, 
b 2 — uj = 0.01, di = —1, d 2 = —1.1, 72 = .1, ipi = vr, ^2 = 0. 



shifts towards the forcing frequency with increasing response amplitude ((bj — u)dj < 0). In this case, the system 
may jump from a small amplitude response to a larger amplitude one with an increase in the forcing. This follows 
from the observation that the instabilities, if they occur, come in pairs and correspond to saddle-node bifurcations 
along the solution branches Rj(F); see Fig. [3Jb) for examples. The necessary and sufficient condition for a pair of 
saddle-node bifurcations to occur is 



a jCj + (6, - u)dj <-\J\{a] + fa - uj) 2 )(cj + d]) , (14) 

for j — 1 and/or j — 2. Note that the prediction that hysteresis can occur for larger detunings is a direct result of the 
dependence of the cell's preferred frequency on response amplitude, an effect which was neglected in previous studies 
for which the imaginary part of the nonlinear coefficient (dj ) was neglected [l(J [H| • 
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FIG. 4: Amplification vs. frequency plots from (J8]l9j) . Amplification was taken to be 72-R1/-F in the single tuning mechanism case 
(a), and R2/F in the double-tuning, uni-directionally coupled cases (b-c). Each curve represents a constant forcing amplitude 
taken from the set F = 10" 3 , 10" 2 5 , 10" 2 Q ...10" S (a) cu = -.01, 61 = 1, Ci = -1, and di = -2, 72 = .01, tp 2 = 0.64tt. (b) 
ai = -.01, a 2 = -.02, 61 = 62 = 1, ci = c 2 = -1, di = -2, d 2 = -1.1, 71 = 0, 72 = 0.01, ip 2 = 0.64tt. (c) ai = -.01, 
a 2 = -.02, 61 = 1, 62 = 1.2, ci = c 2 = -1, di = -2, d 2 = -1.1, 71 = 0, 72 = 0.01, V>2 = 0.64tt. 



IV. RESULTS AND DISCUSSION 



We expect that if each independent amplification mechanism is well-tuned to the forcing frequency Qbj — oj | -C 1) 
and the damping is small the coupled system will have a greater response than a system with only a single tuning 
mechanism. From (| 10[) we see that the leading coefficients ctj for j = 1, .... 7 in © may decrease in magnitude, 
compared to the highest order coefficient ag, as the forcing frequency w approaches the natural frequencies of the 
tuning mechanisms, bj. Figure [4] demonstrates that the amplification, is greatest near the resonant frequency. It 
also makes a comparison between the tuning curves for the coupled system and the system with one amplification 
mechanism suppressed. For the latter case we solve © for with R 2 = 0, and then plot 72^ as a function 

of frequency, for different values of F; the coupling 72 is included so that the same relationship between the hair 
bundle displacement and the magnitude of the transduction current is assumed in both the single and the coupled 
tuning models. Each curve in the diagram shows the variation in response with forcing frequency, holding the signal 
amplitude F constant. Figure HJa) shows that the broadest tuning curve occurs for the loudest sounds. As the 
magnitude of the sound decreases, several effects occur: the amplification increases, the frequency tuning becomes 
sharper, and because dj ^ the preferred frequency shifts as the magnitude of the forcing increases. We used dj < 
for all of the plots, so the preferred frequencies of the oscillation mechanisms decrease as their amplitude increases. 
This phenomenon in turn leads to a region of bistability for driving frequencies u> < bj as described in Section IIIII 
This bistability is the source of the prominent shoulder in the tuning curves that appears for u) < 1. A comparison of 
the height of the peaks in Figure SJ^a) and (b), reveals that amplification of on-resonance forcing is enhanced in the 
coupled system. Also, the frequency tuning curves are sharper and display a smaller shift in frequency with changing 
forcing amplitude in the coupled system. Each of these properties would be a potential advantage of the coupled 
system. If the preferred frequencies of the two mechanisms are not sufficiently well-tuned to each other, then the 
single resonance peak in Fig. [4jb) may split into two peaks as seen in Fig. BJc). 

The algebraic relationship between the magnitude of the forcing and the magnitude of the response of the electrical 
resonance mechanism given in @ allows for experimentally testable predictions to be made. In an idealized situation, 
for which each mechanism is perfectly tuned to the forcing frequency (b\ = b 2 = to) and situated at the Hopf bifurcation 
(ai — a 2 =0), equation §§§ indicates that the response of the electrical resonance mechanism is proportional to 
Fv . The exponent, 5, of this response-versus-forcing relationship, R 2 oc F s , provides a measure of the quality of the 
amplification; small amplitude signals are amplified to a greater extent for smaller values of S. For a system employing 
only a single tu ning mechanism, critically tuned to the bifurcation point and forcing frequency, the expected scaling 
is R oc Fs [nj lll| The smaller exponent of S = 1/9, associated with the well-tuned coupled model, can provide 
more powerful amplification (R 2 /F) than a system with an isolated tuning mechanism. In more realistic situations 
for which some damping and/or detuning is present (dj < 0, bj ^ cj), equation @ predicts a transition from a regime 
in which R 2 oc F for smaller forcings and a regime in which R 2 oc F 1 / 9 for sufficiently large signals. By measuring the 
response of the electrical resonance mechanism for different amplitudes of signal, it is possible to estimate 8, although, 
depending upon which portion of the response-versus-forcing curve is sampled, different estimates for 5 might be 
obtained experimentally. 

In comparing the scaling predictions of our model to experiment, we turn to auditory nerve data. Changes in the 
membrane potential of the hair cell body result in the release of neurotransmitters at the hair cell-auditory nerve 
synapse. Larger depolarizations result in larger amounts of neurotransmitter release and thus a faster firing rate in 
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the auditory nerve. There is biological evidence that, even for nonresonant stimuli, the firing rate at the auditory 
nerve is a nonlinear function of the sound stimulus. However, this nonlincarity, which occurs either due to synaptic 
effects or the relationship between hair bundle displacement and DC receptor potential, is factored out during the 
data analysis, allowing an estimate of S to be made (see discussion in [25l [26j). 

Numerous experiments have been performed in order to estimate the exponent, 6, of the response-versus-forcing 
relationship from experimental auditory nerve recordings. Table I summarizes some recent measurements of 5 taken 
from experiments in both mammals and nonmammals. In each case, S is commonly measured to be smaller than 1/3. 
In mammals, as in nonmammals, two different mechanisms have been proposed to explain amplification, the first due 
to active motion of outer hair cells [U HI] , and the second due to the active motion of the hair bundles [g, H^, [3(| HH . 
While our results do not directly apply to mammalian data, as the coupling of the two active elements may be more 
complicated (32j . it is noteworthy that compression estimates in mammals are similar to those in nonmammals. 

Table I. Experimental estimates of 5 
Nonmammals : 

Owls: between 0.05 and 0.55, with the majority of data points lying between 0.1 and 0.3 [33] 

Pigeons: between 0.22 and 0.6 [34] 

Mammals: 

Guinea-pigs: between 0.2 and 0.25 [25]. 

Guinea-pigs: approximately 0.6 for two low- frequency (1.8 and 2.7 kHz.) fibers, 

and approximately 0.1 for medium- (5.5-6.3 kHz) and high-frequency (20.5-23 kHz) fibers. 

For fibers tuned above 4 kHz, the mean exponent was 0.13 with a standard deviation of 0.04 [35]. 

Chinchilla: direct basilar membrane measurements yield 8 values between 0.2 and 0.7 [36]. 

The observation of S values smaller than 1/3 is interesting in that it indicates that the auditory system achieves 
greater compressions than would be expected from a system with a single tuning mechanism associated with a Hopf 
bifurcation. Experimental measurements of the forcing-versus-displacement relationship for individual hair bundles 
satisfy the R cx F 1 / 3 scaling law expected for a system tuned near a generic Hopf bifurcation [13, EH- Because effects 
at the synapse are removed during the data analysis, any additional compression must occur due to the interaction 
with the electrical resonance mechanism. 

Some models have explained this increased compression by assuming within their model, that the leading nonlinear 
terms are higher than cubic (e.g. [37]). We propose that a more physically motivated way of achieving higher order 
compression would be through the coupling of two systems tuned near a Hopf bifurcation. With the exception of 
three data points, all <5-values in Koppl and Yates' owl data are greater than .1, with the majority of measurements 
lying between .1 and .3 37]. So their data is not inconsistent with what would be expected from a coupled system, 
which at best produces a 5- value of 1/9. If it is the case that the observed enhanced amplification occurs due to 
coupling between the two tuning mechanisms, this is an interesting result, because our analysis indicates that the 
tuning mechanisms must maintain themselves close to the bifurcation point and close to the same frequency for 
maximum amplification to be observed. It would be interesting both from a biological and mathematical perspective 
to understand how such fine tuning is achieved. Some studies have suggested that stochastic effects may help the 
system adjust itself to the bifurcation point [38, 39]. Another has suggested that with certain assumptions about the 
evolution of the bifurcation parameter, self-tuning occurs automatically [io| . Biologically, it has been suggested that 
adjustments in the tension of the hair bundle due to an actin myosin mechanism may act to keep the hair bundle 
properly tuned [H,[4l|. 

Another prediction of the model is that the response of the electrical resonance mechanism may be hysteretic for 
sufficiently large frequency detuning. It is worth noting that this feature arises as a direct result of shift in the natural 
frequency of the tuning mechanisms with increasing response, a feature that was not taken into account in previous 
models. It would be interesting if such bistability could be observed in experimental auditory nerve data. 
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V. APPENDIX I 
Model Proposed by Hudspeth and Lewis 
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The following seven-dimensional model of a bullfrog saccular hair cell was proposed by Hudspeth and Lewis [12], ll3| ■ 
It is a Hodgkin and Huxley type model with three ion channels included in the model, a voltage-gated calcium ion 
channel, a calcium-gated potassium ion channel, and a passive leakage channel. The passive leakage channel is always 
open. The calcium ion channel opens in response to depolarization. The gating of the potassium ion channel requires 
both the binding of calcium to the interior of the channel and an adequate depolarization. The first equation models 
the rate of change of the membrane potential (V m ) due to currents entering the cell through both the transduction 
channels, the calcium ion channels, and the calcium-gated potassium channels. The other six equations account 
for changes in the concentration of internal calcium (Co) close to the membrane, the fraction of open calcium ion 
channels (to), and the fraction of potassium ion channels in each of three closed configurations (Co,Ci,C2) and two 
open configurations (0 2l 3 ). 

Crn^JT = -gc a m 3 (V m - E Ca ) 

at 

9K{Ca) 

(0 2 +0 3 )(V m -E K ) 
-9L{V m -E L )+I 
dCa _ Ug C am 3 (V m - E Ca ) 



dt zFvceiit; 
^ = p(V m )(l-m)-a(Vm)m 

dC 
dt 
dCi 

~~dT 

dC 2 

dt 
d0 2 
dt 

3 = 1 — Co — C± — C2 — 2 



- K,Ca 



k-iCi - hCaCo (15) 
hCaCo + k- 2 C 2 - (k-i + k 2 Ca)Ci 
k 2 Cad + a C oe- Vm/Vaa 2 - (fc_ 2 + C )C 2 
PcC 2 + k_ 3 3 - (a co e~ Vm/v - + k 3 Ca)0 2 



where 



a(V m ) = a e-^ v ^ v - + K A 
(3{V m ) = Poe (v ™ +v « )/VB + K B 
, fc-i 
Ki e ht 

k 2 = ^ 
K 20 e rt 

, fc-3 

3 _ S 3 zFV m 

K 30 e rt 

The amount of current that is injected into the cell, /, is used as the control parameter. All other parameters 
are set to the value used in Hudspeth and Lewis's original paper: g Ca — 4.14 x 10~ 9 S, E Ca = -1 V, gK Ca = 
16.8 x 10~ 9 S, E K = -.08 V, E L = -.03 V, g L = lO" 9 S, C m = 15 x lO" 12 F, U = .02,z = 2, F = 96485.309 C/mol, 
v ceU = 1.25 x 10- 12 L, £ = 3.4 x 10" 5 , K s = 2800 s" 1 , /3 = .97 s" 1 , V = .07 V, V B = .00617 V, K B = 940 s" 1 , 
K A = 510 s-\ V A = -00801 V, a = 22800 s~\ a C0 = 450 s" 1 , T = 295 K, R = 8.314510 ^rp, 5i = .2, S 2 = 0, 
5 3 = .2, /3c = 1000 s-\ k- x = 300 s" 1 , fc_ 2 = 5000 s" 1 , fc_ 3 = 1500 s" 1 , K w = 6 x 10~ 6 M, K 2a = 45 x 10~ 6 M, 
K 30 = 20 x 10 ~ 6 M, and V aa = .033 V. 



VI. APPENDIX II 



Here we present some details of the reduction of the Hudspeth and Lewis model (fT5|) to the normal form ([3]) . Our 
analysis is valid in a neighborhood of the equilibrium solution (V^, Co*, to*, Cq , C^, C 2 , 2 ) for command currents 



12 



I sufficiently close to the critical current I*. In particular, we find that V* z pa —0.04888 V, Ca* pa 1.623 x 10 5 , 
m* re 0.3115, w 0.08429, C£ pa 0.4919, pa 0.1774, O* Pa 0.08960, and C% w 0.1568, for J = 7* Pa 91.3 x 10" 12 A. 

We first translate the fixed point to the origin and nondimensionalize variables as follows: V m — V^(l + X±), 
Ca = Ca*(l + X 2 ), m = m*(l + X 3 ), C = C* *(l + X 4 ), C x = Cf(l + X 5 ), C 2 = C|(l + A 6 ), 2 = 05(1 + X 7 ). 
Moreover, we let I = J*(l + AJ), where A/ measures a small deviation from the critical command current associated 
with the Hopf bifurcation, and allow for a small periodic forcing through the leakage current conductance by setting 
cjl = .92(1 + A<7£(i)). Here Agi,(t) = Ag^if + T) captures the purely oscillatory part of <?l, while g* L is the mean 
conductance; we set g* L = 10~ 9 S, as in the original Hudspeth and Lewis model. The period T is related to the 
forcing frequency 10 in the usual fashion, T = 2tt/lu. The governing equations, expressed in terms of the dimensionless 
vector- valued variable X = (Xi,X 2 , ■ ■ ■ ,Xj) and the parameters AI and Agi,(t), are written 

^=H(X;A/,A 3L (t)) . (16) 

The Hopf bifurcation occurs at X = 0, AI = for Ag L (t) = 0, so H(0; 0, 0) = and the Jacobian matrix DH[0; 0, 0) 
has a pair of purely imaginary eigenvalues ±iu>o, ~ 938 s~ x , with associated complex eigenvectors U, U. The 
remaining eigenvalues Ai,...,As (where A5 = A4) all have real parts less than —2000 (see Fig. [2|), and are 
associated with eigenvectors Vi, . . . , Vg (where, again, V5 = V4). Our convention is to normalize all eigenvectors to 
1. Finally, we diagonalize the linearized problem at the bifurcation point by letting X(t) = z(£)U + z(£)U + ?/i(t)Vi + 
V2(t)~V2 + • • • 2/s(i)Vs, and write the governing equations in terms of the new variables z,z,y= (y±, • • • , j/5) as follows: 

dz 

— = iuj Q z + N z (z,z,y;AI,Ag L (t)), 
dz 

— = -iu; z + N z (z,z,y;AI,Ag L (t)), (17) 
at 

^ = Ay + N y (z,z,y;A/,A 3L (i)) . 

Here A is a diagonal matrix, with eigenvalues Ai, . . . , A5 on the diagonal, and N z and N y contain the nonlinear terms 
in z,z and y, as well as all terms involving the parameters AI and Ag^(t). 

We now use perturbation theory to derive the normal form of the bifurcation problem. Toward this end we introduce 
a small book-keeping parameter e (|e| -C 1) that is a measure of proximity to the Hopf bifurcation. Specifically, we 
let AI = e 2 /i and seek small amplitude solutions of the form 

z(t) = ezx(t, T) + e 2 z 2 {t, T) + e 3 z 3 (t, T) + • • • , (18) 
y(t) = e 2 y 2 (t,T) + e 3 y 3 (t,T) + ■ ■ ■ , 

where T — e 2 t is a slow time variable that captures the slow decay to the oscillatory solutions associated with a 
Hopf bifurcation. Finally, we make some additional assumptions about the magnitude and frequency of the applied 
periodic forcing by letting Agi,(t) = e 3 f(t) and too = u> + e 2 c£>, i.e. we consider small nearly-resonant periodic forcing. 
Inserting these ansatz in equation (|17[) . and expanding in powers of e, we recover at 0(e) 

dzi . . . 

with solution 

zi = A{T)e luJt . (20) 

Here the complex amplitude A(T) satisfies an equation to be determined at higher order. 
At C(e 2 ), we find 

iojz 2 + a z fi + f3 z zf+j z zf + S z \z 1 \' i , (21) 



Of 

where a z ps —756 + 52i, (3 Z pa —86 + 2i, -f z pa 20 — 20i, 5 Z pa —80 — 64i, as well as similar equations for the fast-time 
evolution of y 2 . The general solution of (|2"Tj) is 

Z2( ^ T) = fa + ^jjP) _ lh^ e ^ + nU^ e - 2 ,. t + H T)e^ . (22) 

LO 3w 
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Here B(T) is arbitrary and we set it to zero in the remainder of our computations, since, without loss of generality, 
it may be absorbed by the 0(e) solution (|20[) . The equations for each of the components of y2, which have a similar 
structure to (f2"Tj) . yield solutions of the same form as (|22|) . although with B{T)e lu>t replaced by rapidly decaying 
solutions Cj(T)e Xjt of their associated homogeneous problems. 
Finally, at C(e 3 ), we find 

^ = iLuz 3 + ( iuA - ^— + (45 + 97i)M - (46 + m)\A\ 2 A - (156 - lliUiV"'* + • • ■ • (23) 
at \ oT J 

Here the ellipsis indicates additional terms proportional to e~ tuJt and e ±3lwt , and /i is the coefficient of the e iujt term 
in the Fourier expansion of the periodic forcing function f(t). We have written explicitly only those terms on the 
right-hand-side that are resonant with the solution of the linear homogeneous problem. In order for our perturbation 
expansion to remain valid, this resonant forcing term, which leads to unbounded growth, must vanish. In this fashion 
we obtain the following evolution equation for A(T): 

dA 

— =iuA + (45 + 97i)nA - (46 + m)\A\ 2 A - (156 - lli)/i . (24) 
oT 

Finally, we re-write the equation in the normal form Q by letting 

A(t) = eA(eH)e iut+iv , (25) 
where the phase tp is specified below. We then find 

dA 

— = (a + ib)A + (c + id)\A\ 2 A + Fe w \ (26) 

where 

a + ib = 45A7 + i(uj Q + 97AI), 
c + id = -(46 + 49*), 

F = -(156- lli)Ag LA e i,p . (27) 

Here Ag L l is the coefficient of the e lut term in the Fourier expansion of A<?£, and we choose tp so that F is real and 
positive, i.e. so that 

F= |(156 -lli)A^ Lil | . (28) 

Next, we will consider how two systems tuned near a Hopf bifurcation would interact. Under the assumption that 
the coupling between the two systems is weak and linear, it's straightforward to show that the coupled system can be 
described by the following set of coupled amplitude equations. 

rj A, 

— - = (oi + b 1 i)A 1 + (c x + d 1 i)\A 1 \ 2 A 1 + 7l e^M 2 + Fe^\ (29) 

^ = (oa + b 2 i)A 2 + (pa + d 2 i)\A 2 \ 2 A 2 + l2 e^A 2 , (30) 
at 

where the coefficients of the coupling terms, 71 and 72, are taken to be of order e 2 indicating weak coupling between 
the two systems. Additionally, for coupling between the two systems to occur, it is necessary to assume that the 
resonance frequency of each system is within e 2 of the forcing frequency uj. As with the previous amplitude equations, 
the resonant forcing amplitude, F , is assumed to be small on the order of e 3 and each system is assumed to be tuned 
sufficiently close to the Hopf bifurcation, requiring that both a\ and a 2 be of order e 2 . 
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